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Numerical calculation of the Landauer conductance through an interacting electron 

system in the Hartree-Fock approximation 

Yoichi Asada 

Department of Physics, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152-8551, Japan 

We develop a new numerical method to calculate the Landauer conductance through an 
interacting electron system in the first order perturbation or in the self-consistent Hartree- 
Fock approximation. It is applied to one and two dimensional systems with nearest-neighbor 
electron-electron interaction. 
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1. Introduction 

The interplay between disorder and electron-electron 
interaction in quantum electron transport phenomena is 
one of the most challenging problems. Although the non- 
interacting approximation successfully explains many as- 
pects of experiments, 1 we may need to take into account 
the electron-electron interaction to understand some phe- 
nomena, such as the metallic behavior in two dimen- 
sional (2D) systems in Si-MOS and heterostructures, 2 
and the critical phenomena of the 3D metal-insulator 
transitions. 3 In theoretical works, important corrections 
due to electron-electron interaction in disordered elec- 
tron systems were found. Furthermore, the study of 
the nonlinear a model suggested new universality classes 
for the metal-insulator transition. 9 

The Hartree-Fock (HF) approximation was employed 
in some numerical works to study interacting disordered 
electron systems. 10-18 It enables us to simulate relatively 
large systems. Even at the level of the HF approxima- 
tion, however, our understanding is not yet complete. 
One of the reasons is the lack of numerical simulations 
of the Landauer conductance. As used in the scaling 
theory of Anderson localization 19 (See Refs. 20-22 for 
more detailed discussions on the scaling hypothesis of 
the Landauer conductance.), the conductance is one of 
the most important physical quantities to characterize 
a disordered electron system. We expect that numerical 
simulation of the Landauer conductance would improve 
our understanding on the interplay between disorder and 
electron-electron interaction. 

Motivated by this, we have decided to perform a nu- 
merical calculation of the Landauer conductance in in- 
teracting disordered electron systems in the HF approx- 
imation. As a first step toward it we have developed a 
numerical method, which we report here. 

The generalization of the Landauer approach to inter- 
acting electron systems is a very active topic not only 
to study interacting disordered electron systems but also 
to study transport phenomena through low dimensional 
electron systems. Although much progress has been made 
in the generalization of the Landauer approach, 23-36 the 
calculation method is still being improved. 

This paper is organized as follows: In § 2 the model 
considered is described and in § 3 the Landauer formula 



is described. The HF approximation is described in § 4. 
In § 5, we explain a new numerical method, which we 
call wide band method. In § 6 and § 7, we apply the 
wide band method to ID and 2D systems of interacting 
electrons. The last section is devoted to summary and 
discussion. 

2. Model 

We consider spinless electrons on a 2D square lattice. 
As illustrated in Fig. 1, the system considered consists 
of a sample region with size L s x L y (denoted by S) and 
two semi- infinite leads with width L y (denoted by C). 
We impose fixed boundary conditions in the transverse 
direction. We suppose that electrons are interacting in 
the sample region, while non-interacting in the lead re- 
gion. In the present paper we do not consider random 
potential for simplicity. We take the x direction as the 
current direction and y direction as the transverse direc- 
tion. 

The tight binding Hamiltonian is given by 



H 



H s + He + He 



H v 



(D 



where H s is the non-interacting part for the sample re- 
gion iS, He the Hamiltonian for the perfect leads at the 
left and right, He s the coupling between the sample and 
leads, and H M electron-electron interaction in the sample 
region. They are given by 



H s = — t; 



Ha 



E 

(i,j)(i,jeS) 



cjc 



E c \c] 
(i,j)(i,jeC) 

L I, 



(2) 



(3) 



His = 



E { c ky Cl >y + c i-,y c ^y) 



,(4) 

3/=l 

Here c\(ci) denotes the creation (annihilation) operator 
of an electron at the site i, is the interaction between 
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Non-interacting perfect lead Sample region Non-interacting perfect lead 

Fig. 1. An example of the Landauer geometry to calculate the 
two terminal conductance. Two non-interacting perfect leads arc 
attached to the sample. 



electrons at the sites i and j, and K is a positive uniform 
charge. Hopping is restricted to nearest-neighbors. We 
suppose that the hopping parameters t 8 , te, and tg s are 
real, so the system has time reversal symmetry. 

3. Landauer conductance 

The Landauer conductance g at zero temperature in 
the linear response regime is expressed using Green's 
functions as 

^2 



9 = J'" 



r( L V)G? L >)r( R V)G2 sl ( M ) . (6) 



Here /x is the chemical potential of the system, G r L 1 
and G\ L the submatrices of the retarded and advanced 
Green's functions in equilibrium, and r( L,R ) the matrices 
determined by the attached leads. 23 The matrices ]?( L ' R ) 
are explicitly shown in Appendix . 

It has been shown in Rcfs. 29 and 31 that the expres- 
sion (6), which is well known for non-interacting elec- 
trons, 23 is valid even when electrons are interacting in 
the sample region at zero temperature. Within the HF 
approximation the expression (6) can also be justified 
based on the Keldysh Green's function method. 23 

In the following sections, we only consider the retarded 
Green's function G r since the advanced Green's function 
G a is simply Hermitian conjugate of the retarded Green's 
function. 

4. Hartree-Fock approximation 

We describe the HF approximation in the Green's 
function formalism 37 for the system (1). 

The Dyson equation for G r in the HF approximation is 
obtained by expanding the time ordered Green's function 
in H u and by performing an analytic continuation. We 
have 



G r {e) = G r (e) + G^e)^ 11 ^ G r (e) 



(7) 



Here Gq is the retarded Green's function for the non- 
interacting part of the Hamiltonian Ho = H s +He + He s , 
and 5]( HF ) is the self-energy due to electron-electron in- 
teraction in the HF approximation, which is non-zero 
only in the sample region. It consists of Hartree and ex- 
change contributions: 



.(HF) 



: £ U, 

r;(HF) 



1 f 

— delmG^(e) -K 

J — oo 



delmGUe), 



(8) 



(9) 



where i,j G S. (In more general expression, ImG r (e) 
is replaced with (-i/2) [G r (e) - G a (e)}. For our model, 
they are the same.) From (7), we have 



G r (e) 



H - £( HF ) + 



IT) 



(10) 



Here Hq is the single particle Hamiltonian in the matrix 
form corresponding to the non-interacting part Ho and 
7] is an infinitesimal positive number. The HF self-energy 
E( HF ) is a solution of the self-consistent equations (8)- 
(10). 

In practice there are two difficulties we need to solve: 

• Since the Landauer geometry corresponds to an 
open system, the size of the matrix H is infinite. 
We need to make the matrix size finite to perform 
numerical simulations. 

• To calculate the HF self-energy £( HF ), we need to 
perform an integral over e. 

The first problem can be solved by taking account of 
the effects of the semi-infinite leads in terms of a self- 
energy. 23,3 For example, if we expand the Green's 
function in Ht s in addition to H Ul we obtain the Dyson 
equation for the retarded Green's function G\j G S) 
in the HF approximation in a finite size matrix form that 
is closed in the sample region S. The Green's function 
G\ ■, with i,j G S, is written as 



G r {e) =\e-H s - S^ r (e) - ^ (liF) 



177 



(11) 



Here H s is the single particle Hamiltonian in the matrix 
form corresponding to H s , £( HF ) is the HF self-energy 
which is given by Eqs. (8) and (9), and Y,^ r is the re- 
tarded self-energy due to the attached semi-infinite leads. 
An clement of the self-energy S^^,^, is non-zero only 
when x = x' = 1 or x = x' — L s . The non-zero elements 
are written as (see Appendix ) 



^Xy,Xy'\ e ) 



j£x«(y)xn(y')c(^-A n ), (12) 



where X = 1 or L s . The Green's function G£ • (i, j G S) 
of (11) is exactly the same as that of (10) since we have 
taken into account all orders in Hi s . Now in (11) the size 
of the matrices is finite, L s L y x L s L y , so it is possible to 
perform numerical calculations in principle. 

As for the second problem, a numerical integration was 
employed in similar approaches. 39. The numerical in- 
tegration of the Green's function is not very difficult if a 
simple system is considered. However the numerical inte- 
gration is troublesome in general, so we have developed a 
method to avoid it. The method is explained in the next 
section. 

5. Wide band method 

We have developed a new method, which we call wide 
band method, to avoid the numerical integration as fol- 
lows. 

We change the system considered from the original 
Landauer system, illustrated in Fig. 1 , to another system 
as illustrated in Fig. 2. We call it wide band system. The 
lead region is divided into two parts: the region C up to 
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Fig. 2. The system envisaged in the wide band method. Each 
non-interacting lead is divided into two regions C and W. 



a length L c on both sides of the sample, and the region 
W consisting of two semi-infinite regions. The wide band 
system is described by the following Hamiltonian, 



Tt = H s + He + His + H u . 



(13) 



Here H s , He Sl and H u are already given in (2), (4), and 
(5). The other term He for the leads is made of three 
terms 



He =H C + H w + H c 
which are given by 

H c = -te 
H w 



C i C j' 



e w c \ c i 

i(»ew) 



E 



(14) 



(15) 



(16) 



H c 



tew ^ ^c5 f - L _ 1 j / Cx L ,j/ + c X L ,j/ C *L-l,3/ 
V 

+ C X K ,y C X R +l,y + C XK+1>y C XK ,y^ ■ (17) 

with Xl = —L c + 1 and Ar = L s + L c . Here H c is the 
Hamiltonian for the finite regions C on the sample, H w 
the Hamiltonian for the semi-infinite regions W, and H cw 
the coupling between the regions C and VV. Two hopping 
parameters t w , t cw and one parameter for uniform po- 
tential e w are introduced for the wide band system. If 
t w = tcw = te and e w = 0, the Hamiltonian (13) is ex- 
actly the same as the original Hamiltonian (1), which 
is the system we want to solve. However, as we explain 
below, we take a limit t w ,t cw ,e w — ► oo under certain 
conditions (conditions (24) and (25)). We call this limit 
wide band limit since the band width in the region W 
becomes infinity in the limit t w — ► oo. 

We expand the Green's functions in H cw and H u - The 
retarded Green's function for the wide band system (13) 
in the HF approximation is written as, 



G r {e) =\e-H t - t^ r {e) - £ (HF) + i 



n -1 



IT] 



(18) 



Here H t is the Hamiltonian in a matrix form correspond- 
ing to H t =H s + H ts + H c , S (HF) the self-energy in the 
HF approximation for the system (13) 

1 '' delmG^.(e) 



E„ (HF) 



E * 



K 



(19) 



.(HF) 



7T 



deImG^,(e), 



(20) 



where i, j G S, and E^) 1 " the retarded self-energy due to 
the semi-infinite region W. An clement ^^L, , is non- 

xy : x y 

zero only when x = x = Xl or x = x = Xr . The 
non-zero elements are given by 



"Xy,Xy 



i 2 
tw 



Xn(y)xn(y'K 



e - e y 



-An , (21) 



where X = X^ or Xr. The size of matrices in (18) is 
finite, (L s + 2L c )L y x (L s + 2L c )L y . We use the tilde to 
denote that they are the Green's function and the self- 
energy not for the original system but for the wide band 
system. 

In the semi-infinite regions W we take the wide band 
limit. In this limit the self-energy T,^ 7 " becomes indepen- 
dent of e and the non-zero elements of & w '> r are equal 
to 

^XW = ** E Xn(y)xn(y')( (f - A„) . (22) 

This limit is obtained by taking the limit, 



tw>tcwi^-w * OO) 



while keeping 



t C w/tw — 



v cwl v w 
^w /tw 



= -fi/tt 



(23) 

(24) 
(25) 



(When n = 0, we do not need to introduce e w .) 

A similar idea to use an energy independent self-energy 
can be seen in many papers, for example, Refs. 31 and 
41. Two important differences from previous works are: 

• We keep non-interacting regions up to a length L c 
on both sides of the sample to reduce artifacts of 
taking the wide band limit. 

• The self-energy for the wide band region is chosen 
so that electrons at e = fi are not scattered at the 
boundaries between C and W. This makes it easier 
to reduce the artifacts. 

The self-energy (22) in the wide band limit is equal to 
the self-energy at e = /i for the original system, i.e., the 
self-energy (21) with t w = t cw = te and e w = 0. This 
means that electrons at the Fermi energy e = fj, are not 
scattered at the boundaries between the regions C and 
W. When electron-electron interaction is neglected, the 
Landauer conductance for the wide band system is ex- 
actly the same as that for the original Landauer system 
since the self-energy only at e = \i is relevant for the 
conductance in non-interacting systems. When we take 
account of the electron-electron interaction, the conduc- 
tances for the wide band system and for the original sys- 
tem are no longer the same, because electrons below the 
Fermi energy affect the motion of electrons at the Fermi 
energy through the HF self-energy. To reduce such ar- 
tifacts of taking the wide band limit on the calculated 
conductance, we keep non-interacting regions up to a 
length L c on both sides of the sample. We expect that 
the artifacts of taking the wide band limit decrease as 
L c increases, and they are finally removed in the limit 
L r — > oo. 
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The absence of boundary scattering for electrons at 
e = /i is important for the efficiency in removing the 
artifacts. If electrons near e = fi are scattered strongly at 
the two boundaries, resonant states due to the boundary 
scattering can be formed near e = /i. In this case, the 
Green's function G r (/x) at the Fermi energy shows larger 
fluctuation as a function of L c , and we need to simulate 
systems with longer L c to remove the artifacts. On the 
other hand, our choice minimize the boundary scattering 
of electrons near the Fermi energy, that makes it easier 
to reduce the artifacts of taking the wide band limit. 

By taking the wide band limit, the self-energy E^) r 
becomes independent of e. We define an effective Hamil- 
tonian by 



Then the Green's function (18) is written as 



e _ #(eff)r + 



17/ 



(26) 



(27) 



Note that H^ cS ^ r is not a Hermitian matrix but is a com- 
plex symmetric matrix since H t and X]( HF ) are real sym- 
metric matrices and E( t0 ) r is a complex symmetric ma- 
trix. The effective Hamiltonian has right and left eigen- 
vectors 



i? (off)r |r„) = q n \r n ), 
(l n \H^ r = (l n \q n . 



(28) 
(29) 



Here q n is an eigenvalue, which is complex in general. 
Since the effective Hamiltonian is a complex symmetric 
matrix, the transpose of the corresponding right eigen- 
vector is the left eigenvector, 



(30) 



For convenience, we impose the following normalization 
conditions. 

iln\^rn) — &n.m- (31) 

Then they satisfies the completeness relation, 

X>n><in| = l. ( 32 ) 

n 

By using eigenvalues q n and right eigenvectors \r n ), the 
retarded Green's function is expressed as 23 ' 42 



GUe) = 



y 

^ £ - On + ib n + 17] 



(33) 



Here a n and (—b n ) are the real part and the imaginary 
part of the eigenvalue, q n = a n -ib n , and (f> n (j) = (j\r n ). 

To calculate the self-energy £( HF ), we need to perform 
the integral of the imaginary part of the retarded Green's 
function 

,.( M ,-e c ) = - f delmG^.(e). (34) 

J — e n 



Ji,j ( 



For a moment, we introduce a cutoff parameter e c . We 
will take the limit e c — > oo later. By using the expression 



(33), the integral (34) can be performed analytically as 
Ji,j(fJ; -e c ) = ^2 ) Re [°n{p) - n (-e c )] 

n { 

COS0 n ({l) 



+Im [</> n (i)<j) n {j)]]n 



cos 9 n (-e c ) 



where 9 n (e) (g [— 7r/2 7r/2]) is defined by 



n {e) =Tan" 



q 



(35) 



(36) 



For the second term in (35), we should not take the limit 
e c — > oo before taking the summation over n because of 
the logarithmic divergence 

ln[cos0„(-e c )] = - lne c + In(6„ + rj) + Oie- 1 ). (37) 

The logarithmic divergence disappears when we take the 
summation over n because (32) implies 



n 

Therefore, in the limit e c — > oo we have 
J itj (n, -oo) = Re [<j) n {i)(f>n{j)} n (n) + ^ 



(38) 



+Im [<j) n (i)(j> n (j)]\n 



cos Q n {n) 



b n +V 



>. (39) 



Thus we calculate the self-energy from the eigenvalues 
and the right eigenvectors of if ( off ) without performing 
numerical integration. 

Finally the numerical implementation of the wide band 
method is summarized. First we prepare an initial ma- 
trix for £( HF ). Then we calculate the eigenvalues and 
eigenvectors by diagonalizing the effective Hamiltonian 
H^ eS ' r . From the eigenvalues and right eigenvectors, we 
calculate the self-energy S( HF ), (19) and (20), by us- 
ing (39). We continue this self-consistent iteration un- 
til f]( HF ) converges with enough precision. (To ensure 
convergence, we have used "the method of potential mix- 
ing". 17 ) After convergence, we calculate the conductance 
g(L c ) for the wide band system by using the expression 
(6), where G r ' s (/i) is substituted for G r ' s ([i). We expect 
that the conductance g(L c ) for the wide band system ap- 
proaches to the conductance g for the original system in 
the limit L c — > oo, 



g = lim g(L c ). 



(40) 



In the case of ID system (§6) we extrapolate g(L c ) by 
using an empirical fitting function. In the case of 2D sys- 
tem (§7) we make L c large enough so that the artifacts 
of the wide band limit can be negligible to a good ap- 
proximation. 
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Fig. 3. The conductance g(L c ) for the wide band system calcu- 
lated with the wide band method in ID in the first order pertur- 
bation is shown as a function of L c . We set K = 0.5, (j, = 0.0, 
and U = 0.5. The dotted lines are the fit of (45) to the data with 
L c = [11,30]. The fit deviates from the numerical data when 
L c is much smaller than this range. The solid line indicates the 
conductance for the original system, g 0.9785e 2 //i, obtained 
from (44). 



6. Application to a ID system 

6. 1 System 

Here we calculate the Landauer conductance in a ID 
system, i.e. L y = 1, described by the Hamiltonian (1). 
We suppose that the hopping parameter is uniform in 
the system and set it unity as a unit of energy, t s — 
tt = tg s = 1. We also suppose that two electrons arc 
interacting only when they occupy the nearest neighbor 
sites, 



Ux.x' — 



U (if x' = x ± 1) 
(otherwise). 



(41) 



Transport phenomena in this model was previously stud- 
ied in Refs. 34-36, 43. We set the positive background 
charge K = 0.5 and the Fermi energy /i = 0, which cor- 
responds to half-filling. Since we use the first order per- 
turbation or the self-consistent HF approximation the 
electron-electron interaction should be weak, so we set 
[7 = 0.5. 

6. 2 In the first order perturbation 

First we calculate the Landauer conductance in the 
first order perturbation to test the wide band method. 
Since the system studied here is very simple, the con- 
ductance in the first order perturbation for the original 
Landauer system (1) can be calculated analytically with- 
out making use of the wide band method. By comparing 
the numerical results of the wide band method with the 
analytical results, we test the wide band method. 

The self-energy in the first order perturbation for the 
original system (1) is obtained by replacing G r in (8) and 
(9) with the Green's function Gq for the non-interacting 



Hamiltonian Hq. The first order Hartree term is zero be- 
cause the uniform negative charge of the electrons cancels 
with the uniform positive charge in the background. The 
first order exchange term at \i — is given by 



v (lst) 



- (l<z<L s -l). (42) 



From this self-energy, we obtain the conductance in the 
first order perturbation. When L s is odd, we have a per- 
fect transmission, 



a = 



h 



(L s : odd). 



For even L„ we have 



h 



2(1 + U/ir) 



(L s : even). 



(43) 



(44) 



_{l + U/n) 2 + l_ 

When U = 0.5 and L s is even, we find g ~ 0.9785e 2 //i. 

We have also calculated the conductance in the first 
order perturbation by using the wide band method. The 
numerical result in the first order perturbation has been 
obtained by stopping the self-consistent iteration after 
just one iteration. When L s is odd, the conductance 
is always e 2 /h independent of L c and L s . When L s is 
even, the conductance is reduced from e 2 /h. As shown 
in Fig. 3, g(L c ) oscillates as a function of L c when L s is 
even. 

To remove the artifacts of taking the wide band limit, 
we need to extrapolate the conductance g(L c ) to L c — > 
oo. The extrapolation has been done by using an empir- 
ical fitting function of the form, 

cos(7rL c ) 



g(L c ) = g + a- 



Lt 



(45) 



Here g, a, and y are fitting parameters. We expect that 
the asymptotic value g is equal to the conductance for 
the original system. By fitting numerical data with L c = 
[11, 30], we have found g « 0.9785e 2 //i for any even L s in 
the range L s = [2, 16]. The estimates of a and y weakly 
depend on the range of L c used for the fit. On the other 
hand, the estimate of the asymptotic value g is stable 
against the change of the range of L c . 

The asymptotic value g ~ 0.9785e 2 /ft, for even L s esti- 
mated with the wide band method is in good agreement 
with the analytical result (44) for the original Landauer 
geometry. The difference of the conductance is of order 
10~ 8 to 10~ 6 in units of e 2 /h. The good agreement in- 
dicates that the wide band method works well to esti- 
mate the conductance g by extrapolating the conduc- 
tance g(L c ) for the wide band system to L c — > oo. 

6.3 In the self- consistent Hartree-Fock approximation 

We have then calculated the conductance in ID in the 
HF approximation with the wide band method. When 
L s is odd, we have found g(L c ) = e 2 /h for any L c . When 
L s is even, the conductance is reduced from e 2 /h. We 
extrapolated the conductance for each even L s to L c — ► 
oo by using the fitting function (45), as in the case of the 
first order perturbation. Some numerical data and the 
corresponding fits are shown in Fig. 4. From the fit, we 
have obtained the asymptotic value g for each even L s . 
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Fig. 4. The conductance g(L c ) for the wide band system calcu- 
lated with the wide band method in ID in the HF approxima- 
tion is shown as a function of L c . We set K = 0.5, fi = 0.0, and 
U = 0.5. The dotted lines are the fit of (45) to the data with 
L c = [11, 30]. The fit deviates from the numerical data when L c 
is much smaller than this range. The solid line indicates the con- 
ductance for the original system with L s = 2, g ss 0.9749e 2 /h, 
obtained from (46) and (47). There is no solid line for L s = 16 
since it is too difficult for us to calculate the conductance for the 
original system with L s = 16 in the HF approximation. 



Even in the HF approximation it is not impossible to 
calculate the conductance g for the original system with- 
out using the wide band method. So far we have calcu- 
lated g only for L s = 2. The Hartree term si^P is zero 
when /j, = because of the particle-hole symmetry. The 
exchange self-energy E^ 1 ^ for L s = 2 is a solution of 44 



l-v = -U 



v 2 -l 
2v 2 



1 - -Tan" 



2v 



-1 



+ 



TTV 



(46) 



with v = 1 — sf^. From this equation we find ~ 
—0.1733 for U = 0.5. By using a formula for the conduc- 
tance 



■ 9 = x 



2(1 



.(HF) 
J l,2 



(1 



1 



(47) 



wc find g w 0.9749e 2 //i for L s = 2 and U = 0.5. This 
value of the conductance is indicated with a solid line in 
Fig. 4. The difference between this value and that esti- 
mated with the wide band method is of order 10 _8 e 2 //i, 
indicating that the wide band method works well. 

6.4 Length dependence of the conductance - comparison 
of numerical results 
Figure 5 shows the L s dependence of the conductance 
within the first order perturbation and within the HF ap- 
proximation. For a reference, the conductance calculated 
with the embedding method by R. A. Molina and J.-L. 
Pichard 45 is also shown. In the embedding method, the 
electron-electron interaction is treated exactly, hence it 
is thought that the calculated conductance is also exact. 



Fig. 5. L s dependent oscillation of the conductance g in ID. 
The conductance g in the first order perturbation and in the 
HF approximation are shown as well as that of the embedding 
method. 43 We set K = 0.5, fi = 0.0, and U = 0.5. The data cal- 
culated by the embedding method were provided by R. A. Molina 
and J.-L. Pichard. 45 The lines are a guide to the eye only. 



In all three cases, the perfect conductance g = e 2 /h 
is obtained when L s is odd, and the conductance is re- 
duced from it when L s is even. This even-odd oscillation 
was found in Ref. 43. A similar parity oscillation in the 
Hubbard chain was reported in Ref. 30. 

Within the first order perturbation, the conductance is 
independent of L s for all even L s . However, the result of 
the embedding method, which is thought to be exact, in- 
dicates that the conductance for even L s decreases when 
L s increases. The conductance for even L s in the HF ap- 
proximation also shows that the conductance decreases 
with increasing L s . So we find a qualitative agreement 
between the behavior of the conductance in the HF ap- 
proximation and that of the embedding method. Further- 
more, the results indicate that the HF approximation is 
more accurate than the first order perturbation quanti- 
tatively. 

7. Application to a 2D system 

7. 1 System 

The application of the wide band method is not re- 
stricted to ID systems. Here we apply it to a 2D system, 
i.e., a square system (L s = L y ). The hopping parameter 
is supposed to be uniform and set it unity as a unit of en- 
ergy, t s = tt = ti s = 1. We consider the nearest neighbor 
electron-electron interaction in the sample region, 



U, 



1,3 



U (if (i, j) is a n.n pair) 
(otherwise) . 



(48) 



In our simulation, we set the strength of electron-electron 
interaction U — 0.5. We set the chemical potential (i = 
and the positive charge K = 0.5, which corresponds to 
half-filling. 
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Fig. 6. The conductance g(L c ) for the wide band system calcu- 
lated with the wide band method in 2D in the first order pertur- 
bation for (a)L s = 8 and (b)L s = 16 is shown as a function of 
L c . We set K = 0.5, fi = 0.0, and U = 0.5. The dotted lines arc 
a guide to the eye only. The solid lines indicates the conductance 
g for the original system obtained from (49) and (50) without 
using the wide band method. 



7. 2 In the first order perturbation 

First we have calculated the conductance within the 
first order perturbation to test the wide band method in 
2D as it has been tested in ID in §6. 

Within the first order perturbation, we can calculate 
the conductance without using the wide band method 
since the self-energy £( lst ) in the first order perturba- 
tion for the original Landauer system is obtained as fol- 
lows. The Hartree contribution is zero since the positive 
background charge compensate the negative charge of 
electrons. The exchange contribution is given by 



i(lst) 
- J x,y;x+l,y 



2U 



(Ly + 1)tt 



n=l 



;infc(")sin 2 (fc(">y), (49) 



-i(lst) 
■'x,y;x,y+l 



2U 



£ fcW sin (fcWj 



x sin 



k y n) (y + i) 



(50) 



where k y n) = nn/(L s + 1). Wc have calculated S< lst ) 
by taking the summation over n in (49) and (50) nu- 
merically, and then have calculated the corresponding 
conductance g. 

We have also calculated the conductance making use 
of the wide band method in the first order perturbation. 
Figure 6 shows the L c dependence of the conductance 
g(L c ) for the wide band systems with sample size L s = 8 
and L s — 16. For a reference, the conductance g cal- 
culated from (49) and (50) is also shown with the solid 



Fig. 7. The conductance g(L c ) for the wide band system calcu- 
lated with the wide band method in 2D in the HF approximation 
for (a)L 3 = 8 and (b)L s = 16 is shown as a function of L c . We 
set K = 0.5, /i = 0.0, and U = 0.5. The dotted lines are a guide 
to the eye only. 



line. The figure indicates that g(L c ) fluctuates around 
the value g and the fluctuation becomes smaller when 
L c increases. It is reasonable to assume that the influ- 
ence of the wide band limit is, to a good approximation, 
negligible when L c = 30, 



g wp(L c = 30). 



(51) 



We have found that the differences between the value 
g obtained from (49) and (50) and g(L c = 30) are of 
order 10" 5 to 10~ 3 in units of e 2 /h for L s = [2, 16]. This 
indicates that the wide band method also works well in 
2D. 

7.3 In the self- consistent Hartree-Fock approximation 
We have then calculated the conductance in 2D with 

the wide band method in the HF approximation. Figure 7 
shows the L c dependence of the conductance g(L c ) for 
the wide band systems with sample size L s = 8 and L s = 
16. The L c dependence of g(L c ) is qualitatively similar 
to that in the first order perturbation. We again assume 
that the influence of the wide band limit is negligible 
when L c = 30. 

In the range L c = [21,30], the fluctuation of g(L c ) is 
of order 10~ 3 in units of e 2 /h. Assuming that g(L c ) is 
fluctuating around the asymptotic value g, we can con- 
sider the amplitude of the fluctuation as a precision of g 
in the wide band method. 

7.4 Ls dependence of the conductance 

Figure 8 shows L s dependence of the conductance in 
the first order perturbation and in the HF approximation 
in 2D. For a reference, the value N c e 2 /h corresponding 
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Fig. 8. L s dependence of the conductance g in 2D in the first- 
order perturbation and in the HF approximations. We set K = 
0.5, fj, = 0.0, and U = 0.5. For a reference, the conductance 
corresponding to the perfect transmission N c e 2 /h is also shown. 
(At fi = 0, N c = L a in the 2D strip imposed fixed boundary 
conditions in the transverse direction.) The dotted and dashed 
lines are a guide to the eye only. 



to a perfect transmission, with N c being the number of 
propagating channels, is also shown. The conductance in 
the HF approximation tends to be smaller than that in 
the first order perturbation as it is in ID 

8. Summary and Discussion 

We have developed a new numerical method to cal- 
culate the Landaucr conductance through an interacting 
electron system at zero temperature in the first order 
perturbation or in the self-consistent HF approximation. 
A troublesome numerical integration is avoided by tak- 
ing a wide band limit. We can remove the artifacts of 
taking the wide band limit by increasing the length of 
non-interacting region L c kept on both sides of the sam- 
ple. We have applied it to ID and 2D interacting systems. 

The method does not require much CPU time, so it 
permits us to accumulate many samples in studying in- 
teracting disordered systems. Simulation in the presence 
of disorder is left for future. 

The wide band method has an advantage that the di- 
mensionality is not restricted to one. This method can 
be useful when studying various quantum transport phe- 
nomena not only in interacting disordered systems but 
also in quantum dots, quantum point contacts, quantum 
nanowires, atomic chains, and so on. It is possible to con- 
sider spin degree of freedom. It is also possible to take 
account of other effects, such as spin-orbit coupling and 
a magnetic field, by generalizing the wide band method. 
In the presence of such effects, the effective Hamiltonian 
might be no longer a complex symmetric matrix. When 
it is not a complex symmetric matrix, we need to cal- 
culate not only the right eigenvectors but also the left 
eigenvectors. 
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Appendix: Green's functions for semi-infinite 
leads and some related matrices 

In this appendix, the retarded Green's function for 
semi-infinite leads and some related matrices are shown 
explicitly. We suppose that the leads are described by 
the Hamiltonian (3) and coupled with a sample by (4). 

First we calculate the retarded Green's function g r (e) 
for the isolated lead at the left (x < 0). In particular, 
the important elements are those at the surface of the 
lead, i.e., the elements g$ y 0y '( e )- The eigenfunctions and 
eigenvalues of the lead are written as 



where 



Xn(y) 



-sm[k x (x-l)]xn{y), (A-l) 



-It i cosk x +UX n . 



-2cosfcl n 



= sin (k^J 



(A-2) 



(A3) 



(A4) 



Here k x (e [0, n]) is the wave number in the x-direction, 
and k y n ^ — Tm/(L y + 1). Using them we have 

1 ^ v f \ 

ffo»,oi/'(e) = Y t ^2xn(y)xn(y')C - Kj • (A-5) 



Here the function Q(z n 



C(Zn) 



is given by 

{Zn < -2) 



f-H/l-^f (-2<^n<2) (A-6) 



1 



(z n > 2). 



2 V 4 

Similarly, the Green's function for the other lead at the 
right (x > L s + 1) is obtained as 



1 



9 r L e +ly,L e + ly>(e) = l~^2Xn(y)Xn(y')( 



~\n .(A-7) 



A transverse mode n is a propagating mode if \^i/t( — 
A„| < 2 and an evanescent mode otherwise. The Fermi 
wave number k^ (0 < k n < ir) of the propagating mode 
n in the x-direction is determined by 



ik 



n 

4 



(A8) 
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with z n = n/ti - X n . 

When calculating the Landaucr conductance the semi- 
infinite leads and the sample are coupled by (4). In the 
Green's function method the effects of the lead are taken 
account in terms of the self-energy. 23,38 The retarded 



,(*)r 



due to the leads becomes non-zero 



self-energy T, xyx , yl 
at the surfaces of the sample, i.e., when x = x' = 
1 or x — x 1 — L s . The non-zero elements are de- 

tls90yfly>{ e ) an d ^L s y,L s y'( € ) = 

_i y i(e). Using the expression (A-5) we have 

^Zxy'^) = f Y,Xn{y)Xn{y')C - A„) , (A-9) 



fined by E^e) 



ti.g r L, 



where X = 1 or L, 



The L y x L y matrices and r( R ', which ap- 

pear in the Landauer formula (6), are defined by 



i [^l s yx sv >{t) - ^l s yx s y>(t)\- At e = M we have 
2i 2 



t(L,R) 



(m) 



L Y^ Xn{y)Xn(y')sink 



(A-10) 



Here the summation is taken over propagating modes. 
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